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Heat and work for quantum systems governed by dissipative master equations with a time- 
dependent driving field were introduced in the pioneering work of Alicki [J. Phys. A 12, L103 
(1979)]. Alicki's work was in the Schrodinger picture; here we extend these definitions to the 
Heisenberg and interaction pictures. We show that in order to avoid consistency problems, the full 
time derivatives in the definitions for heat flux and power (work flux) should be replaced by partial 
. . . time derivatives. We also present an alternative approach to the partitioning of the energy flux 

' which differs from that of Alicki in that the instantaneous interaction energy with the external field 

' is not included directly. We then proceed to generalize Alicki's definition of power by replacing the 

, original system and its external driving field with a larger, bipartite system, governed by a time- 

independent Ifamiltonian. Using the definition of heat flux and the generalized definition of power, 
^ ' we derive the first law of thermodynamics in differential form, both for the full bipartite system 

Q , and the partially traced subsystems. Although the second law (Clausius formulation) is satisfied 

• for the full bipartite system, we find that in general there is no rigorous formulation of the second 

' law for the partially traced subsystem unless certain additional requirements are met. Once these 

, requirements are satisfied, however, both the Carnot and the Clausius formulations of the second 

law are satisfied. We illustrate this thermodynamic analysis on both the simple Jaynes-Cummings 
model (JCM) and an extended dissipative Jaynes-Cummings model (ED-JCM), which is a model 
for a quantum amplifier. 
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I. INTRODUCTION 



Thermodynamics of quantum systems has intrigued scientists since the early years of the development of quantum 
^^O ' mechanics, von Neumann [l| defined an entropy function which is based on the density matrix, and is similar in 
spirit to the Gibbs entropy. It takes into account both the populations and the coherences of the density matrix. 
^ Born [2] partitioned the energy of various types of quantum statistical systems into heat and work. Tolman Q 
Oh! defined an entropy function which is identical with the Gibbs entropy, and which is based on a quantum version 
^Jj ■ of Boltzmann's if-function. He showed that Carnot's formula holds for a quantum canonical distribution, and that 
Ch Clausius's inequality holds in general. Scovil and Schulz-DuBois [3| analyzed a three-level quantum system coupled 
to two thermal reservoirs. Much later, a master equation for open quantum systems in the weak coupling Markovian 
regime was developed by Lindblad, who showed that under a trace-preserving completely positive map the relative 
entropy is non-increasing [3]. Spohn defined the entropy production function, and showed it is always positive 0]. 
. ^ Pusz and Wornowicz defined work in a general C* algebraic context for systems with varying external forces ^ . Based 
, on the Markovian master equation, Alicki defined work and heat for systems with time dependent Hamiltonians. 
■ The ideas of Alicki were implemented in a series of papers by Kosloff and coworkers on quantum heat engines 

. ; [a, P- 

This work continues the work of Alicki. The main focus of this paper is laying the foundations for a complete 
thermodynamic analysis of bipartite systems, and to apply this novel approach to quantum optical systems. It is 
not obvious how to generalize the thermodynamic definition of Alicki to quantum optical systems, since Alicki's 
definition of work presupposes a time-dependent external field, whereas in quantum optical systems the Hamiltonian 
is time- independent. To solve this problem, we begin by noting that the time-dependence of the Hamiltonian in the 
semiclassical treatment arises from the time- independent Hamiltonian in the quantum system because we are looking 
at only a part of the larger system. For example, in the absence of dissipation the total energy of the bipartite 
matter-light system is conserved, but the energy of each subsystem in general is time-dependent. If one observes the 
energy changing with time in one of the subsystems, there is no way of distinguishing whether this is as a result of 
external forcing or because the subsystem is part of the larger bipartite system. Therefore, if work can be defined 
for the former we can expect that it can be defined for the latter as well. The central result of this paper, eq. 1151 
provides precisely such an expression for the work, expressed in terms of the Hamiltonian of the full bipartite system. 
Inspection of this equation shows that it is not restricted to quantum optical systems, and is in fact completely general, 
applying to any bipartite system. Nevertheless, the application to quantum optical systems is particularly intriguing 
since it lays the foundation for a thermodynamic analysis of light-matter interactions in non-equilibrium systems. 

In Section \n~K[ we review Alicki's definition of quantum heat and work. In Section III Bl we present a minor gen- 
eralization of Alicki's definitions for unipartite systems that provides thermodynamic consistency in the Schrodinger, 
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the Heisenberg, and the interaction pictures. In section Hi CI we suggest an alternative partitioning of energy into heat 
and work that differs from that of Ahcki and presages the bipartite formulation in Section IIIII In Section III Dl we 
analyze the different formulations of the second law for unipartite systems. Section IIIII on bipartite systems, is the 
main part of this paper. We begin with the derivation of the general formula for work flux (power) and heat flux in 
bipartite systems that generalizes our unipartite formulation of the first law in Section [II CI In section fill Bl we show 
that although the second law in the Clausius formulation is satisfied for the full bipartite system, in general there is 
no rigorous formulation of the second law for the partially traced subsystem unless certain additional requirements 
are met. Once these requirements are satisfied, however, both the Carnot and the Clausius formulations of the second 
law are satisfied. In Section HlI CI we apply the general formalism to the simple Jaynes-Cummings model (JCM) and 
to an extended dissipative JCM, which serves as a quantum amplifier. Section [TVl concludes. 



II. UNIPARTITE SYSTEMS 



A. Energy flux of systems governed by master equations 

The master equation for a system coupled to thermal reservoirs in the Schrodinger picture and the Markovian 
regime is given by: 

P^C[p]^Ch[p]+CM. (1) 

where C[p\ is a general Lindblad super operator, Ch[p\ — -~j^[H{t), p] is the Hamiltonian Liouville super operator, and 
Cd[p] is a dissipative Lindblad super operator. Based on the master equation one can calculate the average energy 
by: 

{E) = TT{pit)Hit)}. (2) 

Alicki partitioned the energy into heat and work by taking the time derivative of equation [21 Heat is defined by: 

t 

|Tr|^i7(0}di', (3) 



where work was defined originally by Pusz and Wornowicz [Sj as: 

W^jTr^pin'J^l^t'. (4) 



A consistency problem arises when one applies these definitions in the interaction picture, and one has to be careful 
when defining a Heisenberg picture for a system governed by dissipative evolution. Moreover, according to eq. [4| a 
system whose interaction with some external degree of freedom is governed by a time independent Hamiltonian does 
not perform any work (as in the case of quantum electrodynamics where field-matter interactions are governed by 
time independent Hamiltonians). These two issues will be addressed in the sections to follow this brief introduction. 



B. Heat and work in various pictures for unipartite systems coupled to a time dependent external field 

1. The Schrodinger picture 
The time evolution of the average value of an operator is given by: 

where in the Schrodinger picture dp^ /dt — dp^ /dt and dA^ /dt — dA^ /dt. Consider a system whose evolution is 
governed by equation [TJ The first term in eq. [5l becomes: 

TtI^A^} = Tr{/:[p^]A^} = -'-Tt{[H^,p^]A^} + Tr{/:,[p^]A^| . (6) 
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Using the definition for quantum heat (eq. [3]) the heat flux in the Schrodinger picture is given by: 

Q = Trj^i/^} = Tr l^i/^} = Tr{U[p']H% (7) 
where in the last equality we have used the fact that TT:{Ch[p^]H^} = 0. The power (work flux) is thus given by: 

Note that in case of only Hamiltonian dynamics there is no heat involved, and if the Hamiltonian is time independent 
there is no work done by the system. 



2. The Heisenberg picture 
Consider a system whose evolution is purely Hamiltonian: 

pS^Cu[p'^]^-'-[H^p% (9) 

Eq. [9] can be formally integrated: 

p'{t) = U{t)p{Q)uHt), (10) 

where U{t) is a unitary operator. If the Hamiltonian is time independent then U{t) — exp~Ti^^* . Otherwise, U{t) is 
obtained by the time ordering procedure p^ . The essence of the Heisenberg representation is that operators 'move' 
in time, and it is defined by [l4| : 

A" = U\t)A^U{t), (11) 
where U{t) is the same unitary operator that appears in eq. 1101 Note that in this case is time independent: 

pH = Lu^[H^^p^]U + U^pSU = {{[H",p"] [i^^,P^]) = (12) 

Note also that if is time independent {t) = H^{t) = H^{Q). The average value of an operator is independent 
of representation: 

{A^) = Tr{p^A^} = i:v{U\t)p^U(t)U\t)A'^U{t)} = Ti{p^ A^} = {A"). (13) 

Consider now a system whose evolution is governed by eq. [1] That is, its evolution is also dissipative. A solution 
in the form of eq. [TO] can not be obtained due to the non- unitary nature of In this case p^ is time dependent: 

pH ^ Lu^H^.p^p + U^pSU = U^Cdip'^p EE £f [p^]. (14) 
n 

The evolution of the density operator in the Heisenberg picture is due to the dissipative part in the master equation, 
while the Hamiltonian evolution is still canceled out. The time dependence of p^ is not an artifact of a definition or 
a transformation. Consider any function of the density matrix, for example purity ((p) = Tr{/5^}). It is obvious that 
under unitary evolution, purity does not change with time (as can be shown by a Taylor expansion of the density 
operator). However, under dissipative dynamics purity may change with time, and this should be the case in any 
physical picture. Therefore, when one calculates the average value of an operator there are two contributions to its 
evolution: the 'moving' density operator, and the 'moving' observable. This can be seen by: 

dt dt \ dt j {' dt j ^ ^ 

The first term in eq. [15] is obtained by substitution of eq. [TJ] 

Tr[^AA^Tr{C^[p"]A"}, (16) 
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while the second term is given by: 



Substituting in eq. [TBI and eq. [TTl shows that the definitions for heat flux and power in the Heisenberg picture 
are identical to those in the Schrodinger picture: 

Tr l^i/^l = Tr{C^[p"]H"} = Tr{Cl[p']H^'}^Q (18) 



dt 



-{'''^}^H^''(f)"}-'{4fn-- 



(19) 



Geva and KoslofF rewrite eq. [5] using the cyclic invariance of the trace. By doing this they obtain an alternative 
Heisenberg picture in which the observable is evolved by the dissipative part of the master equation together with the 
Hamiltonian part [IB], [l6j : 



where: 



dt [ dt J I at J L J J 

Sr*^AS]^ t., J „g 9A^ \ _ ( s 



Tr{p*/:*[A*]} + Tr^p*^| =Tr|p — 



and £* [A] is the full Lindblad super operator as a function of A. However, since still appears in eq. [5T] their 
definitions for heat flux and power [l5| are formally equivalent with Alicki's definitions in the Schrodinger picture: 



Q = Tr{p^/:S[^^]} = Tr{/:rf[/]if^} 



where C*^[H] is the dissipative Lindblad super operator as a function of H. 

3. The interaction picture 

Consider the following Hamiltonian: 

H{t) =Ho + V{t), (23) 

where Hq is usually a simple Hamiltonian (for example the eigenenergy Hamiltonian of a system), and V{t) is the 
coupling to some external degree of freedom or field. The interaction representation is defined by: 

A' = U^.m'^Uoit), (24) 

where Uq is the propagator associated with Hq. The average value of an operator is independent of any such trans- 
formation: 

{A^) = Tr{p^A^} = TT{U^,it)pSUom^om'Uoit)} = Tv{p'A'} = (25) 
Differentiation of (A^) yields: 

diA^) _ WA^})^^(d_l^A^^(^,dA^X ^^^^ 



dt dt I (it I I dt 
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Substituting eq. [24] into the two terms in eq. [26] and using the relation Uq = —jrHoUo yields: 

-K'l^A'l = -i:Tv{p'lA'.V]} + Tt{C'jlp']A'} (27) 

^^'^}->'''[''»'-^'i> + i^{'''(w)')' 

where is the dissipative Lindblad super operator in the interaction picture. Substituting into cq. 1271 and 
eq. [28] yields: 

TTf^^H'y-^TT{p'[H^,V']}+Tr{C',[p']H'}^-^T^^^^ (29) 

Tr{/f^} ^ ^r{/[i.„^V^]KTr{/(f )]4Tr{/[i^,^^^]HTr{/(f )]. (30) 

Adding eq. 1291 and eq. [30]yields dE^ /dt. However, each term differs from the corresponding terms in the Schrodinger 
picture by ^i/ftTr{p'^[i7^, V^]} (see eq. [7]and eq. [5]). Thus heat and work should be redefined to avoid a consistency 
problem when one moves from the Schrodinger picture to the interaction picture. We redefine heat and work as 
follows: 

The definitions for heat and work in eq. [3T] differ from Alicki's definitions for heat and work (eq. [3] and eq. [1] 
respectively) only by replacing the full derivatives with respect to time with the partial derivatives with respect to 
time. This replacement is necessary as it avoids the extra commutator term ^i/h'Tv{p^[HQ ,V^W which arises from 
the full time derivative in eq. 1261 



C. Alternative approach to heat and work for unipartite systems coupled to a time dependent external field 

Consider a unipartite system that is governed by eq. [1] and is driven by a time dependent Hamiltonian of the form 
of eq. 1231 The energy flux of such a system is given by substitution of H{t) into eq. [5] in the Schrodinger picture: 

The average value of the Hamiltonian {{E) = Tv{p^ [t)^) contains contributions from both the system energy and 
the interaction energy. Even if the system reaches a steady state (p = 0), the quantity E is not necessarily zero at 
steady state due to the term Tr Let us consider the energy flux of the system alone ((£^o) = Tr{p'^i7(f (i)}): 

^ ^ ^ {¥^4 ^ Tr{/:.[p"]i^on + ^^{Ck[p']hI}. (33) 

We now define heat flux by: 

Qo = Tr{£rf[p]ffo}, (34) 

and power by: 

Po = Tr{/:,[p]i/o} = -'jTy{p[Ho, V{t)]}. (35) 



These definitions are novel to the best of our knowledge, and have the property of being identical in all physical pictures 
(the expression in the Heisenberg and interaction pictures can be obtained easily through the cyclic invariance of the 
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trace). This partitioning of the energy flux captures the spirit of AHcki's partitioning in the sense that power is defined 
via the Hamiltonian superoperator and heat flux is defined via the dissipative Lindblad superoperator (although the 
present definition for Qo involves only the bare Hamiltonian Hq whereas Alicki's definition for Q involves the full 
Hamiltonian H = Ho + V{t)). Moreover, this partitioning of the energy flux has the same structure as the partitioning 
we will use in section HIA (eq. |45]and eq. I46p where the energy of a partially traced system within a bipartite system 
is considered. Finally, note that the quantity (Eq) does reach a steady state since Hq is time independent. 



D. The second law for unipartite systems 

The ultimate test for partitioning energy should be the fulfillment of the second law. The most general formulation of 
the second law is that of Clausius, which states that the total entropy change of a closed system must be greater than or 
equal to zero. In that respect, an open system together with the reservoirs that are coupled to it is in fact closed. Thus 
the entropy changes of the system and of the reservoirs must fulfill the relation: AStotai — system-'r ASj-eservoirs ^ 0. 
For this relation to hold there is no need to define work. When the system operates in a cycle, the only contribution 
to the total entropy change is the entropy change of the reservoirs {AS system = 0), and one can derive a Carnot 
formulation of the second law, which sets a limit on the engine efficiency. However, in order to have a Carnot 
formulation of the second law, work must be defined. 

To derive the Clausius and Carnot formulations of the non-equilibrium second law of thermodynamics in differential 
form, one replaces the equilibrium thermodynamic quantities of heat and work by thermodynamic currents (or fluxes) 
of heat flux and power. In order to satisfy Clausius's formulation of the second law for driven unipartite systems at 
all times we believe that one has to use Alicki's definition for heat flux. However, at steady state, in order to obtain 
Carnot's formula one must use our definitions (especially for power), as we now explain. 



We begin with Spohn's entropy production function |9| for a system coupled to two reservoirs: 

9S Qh Qc f^f,, 



where ^ is the system's entropy change defined via the von Neumann entropy, and — is the reversible entropy 

current from/to the hot (cold) thermal reservoir. Spohn showed that for a completely positive map (such as the 
Lindblad super operator) [§]: 

cr > 0. (37) 

Eq. [37] is Clausius's general formulation of the second law in differential form, which is valid at all times. 

At steady state ^ = 0, and hence we need to analyze only the thermodynamic currents from/to the reservoirs. 

Note that Qh{c) = "^"''{^dH(c)[p]H} has two contributions, namely Qh(c) = Qh{C')o + Qh(c)v, where Qh{C')v = 
TT{CdH{c)[pW} ■ Writing eq. [33] for the case of two reservoirs yields: 

^0 = Qho + QcQ + Po- (38) 

At steady state {Eq — 0), Q^q — —{Qho + Po'')- If assume that = 0, substitution of eq. [35] into eq. [31] 

yields Carnot's famous formula in differential form: 

'n = —p—<^ ■ (39) 

Qho 

Note that a similar Carnot formulation of the second law does not generally exist using Alicki's formalism, due to the 
fact that generally E ^ 0, for example in the non-resonant semiclassical ED-JCM discussed in [iTj . 
For completeness, we note that if Qc{h)v — 0, eq. [35] is trivially satisfied. 



III. THERMODYNAMICS OF BIPARTITE SYSTEMS WITH TIME INDEPENDENT HAMILTONIANS 



A. Energy fluxes for bipartite systems 



Consider a bipartite system whose evolution is governed by eq. [T] The time independent bipartite Hamiltonian 
operates in a C™ (g) C" Hilbert space, and it is given by: 



H — Ha + Hb + Vabi 



(40) 
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where Ha = Ha <^ Is and Hb = 1a ® Hb- For simplicity we consider direct dissipation through subsystem A 
only. However, the analysis that follows can be generalized easily to the case where both subsystems dissipate energy 
directly to the environment. The generic form of the dissipative Lindblad super operator is given by: 

^d[PAB] ^ ^Tkiakilcrl, (TkPAB] + [PAB(Tl,crk]) + bki[(Tk, ctIpab] + [PAB(Tk, (^1])} = ^C'-J'\pab], (41) 

k k 

where <Tfc(<T^) is the kth lowering (raising) operator of subsystem A in the tensor product space (Tfc((T^) — Cfcyi (cr]^^ ) ® 
I-b), Tfc is the kth decay rate, a^, and bk are prefactors whose ratio is a function of the temperature of the kth reservoir. 

The derivation that follows is general for either the Schrodinger, Heisenberg, or interaction representation. For 
convenience we will work in the Schrodinger representation. Using the definition for heat flux (cq. [21]), the energy 
flux of the full bipartite system is given by: 

Eab = Tr{Cd[pAB]H} = Q = Y,Qk, (42) 

k 

where Qk = TtIC'^^ H} is the heat flux associated with coupling to the kth reservoir. Note that since the Hamiltonian 
is time independent there is no work involved. Thus, to an outside observer looking on the bipartite system as a 
whole the full system is only dissipating. Is this the whole story? No. 

Let us examine the energy flux of subsystem A (the subsystem for which direct dissipation occurs): 

Ea = TtaIpaHa} = Ti{C[pab]Ha} - -^Tt{pab[Ha, Vab]} + TT{Cd[pAB]HA}, (43) 

where Ha is the Hamiltonian of subsystem A without the tensor product with 1b, and pA — Tt^b{pab} is the partial 
density matrix of subsystem A. The first equality in equation 1431 stems from the following algebra: 

Ea = TvaIpaHa} = J] pM = E (E Pa^'^H^a 

= E Pa^^^H^aC - E PAh'^Hi^''" = T^aMpabHa} = Tr{C[pAB]HA}. (44) 
Comparing eq. l43land eq. [44] we associate the first term on the RHS of eq. [43] with power: 

Pa = -^Tr{pAB[HA, Vab]}. (45) 
while we associate the second term with heat flux: 

Qa = Tv{Cd[pAB]HA}. (46) 

Before we continue we wish to emphasize the differences between our novel definitions for heat flux and power (eq. 
[35] and eq. [IS] respectively) and these of Alicki (eq. [3] and eq. [J] respectively) . Firstly, our definitions refer to a 
single degree of freedom within a bipartite system. Secondly, the dimensionality of the operators in our definitions is 
m ® n (due to the expansion in eq. I44p as opposed to m when a single degree of freedom is coupled to an external 
field (Alicki's work). Thirdly, power emerges simply from the commutator between Ha and Vabi and not from a time 
dependent field. The last difference is striking since it indicates that the Hamiltonian need not depend on time to get 
work. In fact, whenever an external field is replaced by a quantized degree of freedom, the time dependence of the 
Hamiltonian is removed at the expense of increasing the dimensionality of the Hamiltonian. 
Similarly, we can analyze the energy fiux of subsystem B: 

Eb = ^tb{pbHb] - -'tT:t{pab[Hb, Vab]} = Ps- (47) 
n 

There is no contribution to the energy flux of subsystem B from the dissipative part of the Lindblad super operator 
(this is physically expected since the dissipation is through subsystem A only): 

Ti{Cd[pAB]HB} = -y^rfcTr|(afc([cr^,crfcp^fl] + [pabct^, Tfc]) 

k 

+ hk{\cTk,rTlpAB] + [pABfTk,rTl]))HB] (48) 
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To see this consider just the terms Tr{([<T^, ct^pab] + [pabc^, CfeD-f^s}: 

Tr{{[crl,(TkPAB] + [PAB(Tl,rrk])HB} = T^r{pAB{HBcrlcrk - icT^Hscrk + crlcrkHB)} 

= Tr{pAB[{^A ® Hb){(jIa ® 1b)Ka ® Is) - 2{crl^ (g) ® HB){(JkA ® Is) 

+ (o-t^ ® lB)(crfcA «) 1b)(1a iJs)} = (2 - 2)TT{pABCrlA<^A ® = 0. 

Similarly, Tr{{[(Tk,crlpAB] + [PABcrk, (tI])Hb} = 0. 

The first law of thermodynamics in differential form can be formulated in two different ways. The first way relies 
on looking at each subsystem independently; for subsystem A it is obtained by combining eq. 32] with the definitions 
in eq. SHI and eq. SS] 

Ea^Qa + Pa- (49) 

The first law for subsystem B is given by eq. 1471 

The second way relies on looking at the full bipartite system; it is obtained by substituting eq. [46] for the heat flux 
into eq. |42l 

Eab = Qab = Qa + Qv, (50) 

where Qv = T^'^{^d[PAB]VAB} and Qa was defined in eq. HH] By adding together eq. I^and eq. [IT] we see that in 
general: 

Ea + Eb=Pa + Qa + Pb^Qa + Qv = Eab- (51) 



B. The second law for bipartite systems 



The second law of thermodynamics in differential form for the full bipartite system is once more obtained via 
Spohn's entropy production function defined in [lo| : 

a^^ + J>0, (52) 

where is the entropy production associated with the full bipartite density matrix via differentiation of the von 

Neumann entropy [T] {Ssys — —fcsTrjpyiB In Pab}), and J is the entropy production associated with the reservoirs 
(via the heat flux from/to the reservoirs) and it is given by: 

J=-5lAQfc. (53) 
fe 

Here —PkQk is the differential change in the entropy of the fcth reservoir, and f3k — ^/ksTk where Tk is the temperature 
of the fcth reservoir. Eq. [52] represents the differential form of the second law of thermodynamics in Clausius's 
formulation since the sum of the entropy changes of the system and reservoirs is guaranteed to be positive. 

Is it possible to formulate a second law of thermodynamics for a subsystem within a bipartite system? Consider an 
entropy production function that is based on thermodynamic currents of subsystem A only, which is coupled to two 
thermal reservoirs: 

a c 

aA = ^ + Ja, (54) 

where is the entropy production associated with the partial density matrix via differentiation of the partial von 
Neumann entropy [1,] , and J a is the entropy production associated with the direct coupling between subsystem A and 
the reservoirs: 

Ja = -PhQha - PcQcA, (55) 

where Qha = Tr{£rf/f and similarly for Qca- Note that each heat flux component in J has two contributions, 

e.g. Qh = Qha + Qhv , where Qhv = Tr{£rf/f [p] Vab}. Unlike cr, a a is not necessarily positive due to the oscillatory 
nature of the partial entropies under Hamiltonian evolution, for example, the partial entropies in the Jaynes-Cummings 
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model (JCM) [Tl| [T^. However, due to the dissipative nature of the Lindblad superoperator subsystem A is expected 
to reach a steady state, at which point = and Ea = 0. Moreover, at steady state, (1) if Qhv = Qcv = 
(which implies that J — Ja), and (2) if the main source for entropy production is the dissipation of subsystem A 
(J > ^fp-), then 

erf = = J > 0. (56) 

If the above two conditions arc fulfilled, eq. [56] corresponds to Clausius's formulation of the second law for the 
subsystem A, and Carnot's formula follows from it trivially. Note that condition (1) appeared in our earlier discussion 
of Spohn's entropy production formula for unipartite systems (Section IV). Eq. is verified numerically for the 
quantum ED- JCM in a forthcoming publication [l3|. 

It is interesting to consider the differences between a a for the partially traced system and a for the forced unipartite 
system described in Section fll Dl We expect deviations between these two to be significant when the entropy of the 
subsystems oscillate, since the unipartite system entropy can only increase. Thus, at short times we expect that there 
may be significant differences, while at steady state, when the entropy of the partially traced bipartite system becomes 
constant, we expect the differences to be small. This is borne out by numerical studies, where transient collapses and 
revivals are observed in the partially-traced bi par tite treatment at short times, while only minor deviations from the 
unipartite values are observed at steady state . Note that these are the same oscillations that are responsible for 
the breakdown of the Spohn entropy production formula for the partially-traced bipartite system, as discussed in Eq. 

M 



C. Applications 

1. The Jaynes-Cummings model (JCM) 

The Jaynes-Cummings model is the simplest fully quantum model that describes the interaction between light 
and matter [23| • In this model a two- level material system interacts resonantly with a single quantized cavity mode 
according to the following Hamiltonian (here and throughout the whole section we omit the superscript S): 

H = Hm+Hf + Vmf, (57) 

where Hm = hoJaO'e (8> 1/ is the matter Hamiltonian, Hf = hufla a^a is the field Hamiltonian, and Vmf = fiX[a^ (g) 
a + a~ ® a^] is the coupling Hamiltonian (cr"*" and ct~ are the raising and lowering matter operators, respectively). 
The matrix form of the matter operators is given by: 

























(T = 




I) 




I 





The system is governed by the following master equation in the Schrodinger picture: 

Pmf = Ch[Pmf], (58) 

where Ch[pmf] = —j^[H,Pmf] is the Hamiltonian superoperator. The energy flux of the full matter-field system is 
given by: 

Emf = TT{Ch[Pmf]H} = 0. (59) 

Note that since the evolution is purely unitary the total matter- field energy is constant. Let us examine the energy 
fiux of the matter: 

E,n = TTm{PmHm} = - -TT{Pm.f[Hm., Vmf]} = Pm, (60) 

where iJ„j is the Hamiltonian of the matter without the tensor product with 1 j , and Pm is the power associated with 
the atom. The energy flux of the selected cavity mode (field) is given by: 

Ef = Trf{pfHf} = ^^Tr{pmf[Hf, V^f]} = Pf, (61) 
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where Hf is the Hamntonian of the field without the tensor product with t„i, and Pf is the power associated with 
the field. Generally, E^ + Ej = Tr{pmf[Hm + Hf, Vmf]} ^ h^\TT{pmf{uJm - w/)(cr+ (g) a - cr" «) a^)} ^ E^nf = 0. 
However, under perfect atomic-field resonance {[Hm + Hf, Vmf] — 0): 

Pm = -Pf. (62) 

and hence Em + Ef = E„if — 0. 

No entropy is produced during the interaction: 

(7 = 0. (63) 

This can be understood in two ways. Firstly, under pure Hamiltonian evolution the eigenvalues of the bipartite system 
are a constant of the motion [19] and hence ^ = 0. Secondly, since there is no coupling to an external heat reservoir 
(closed system) no heat is produced and hence J = 0. However, partial entropies do change with time and the entropy 
content of the partial systems gives information with respect to entanglement and disorder [l^ . 



2. The extended dissipative Jaynes-Cummings model (ED-JCM) 



Consider a three-level system interacting resonantly with one quantized cavity mode and two thermal photonic 
reservoirs as depicted in Fig. [1] The system is governed by the following master equation in the Schrodinger picture: 




FIG. 1: Three level system interacting with two heat reservoirs (hot and cold) and a quantized cavity mode. 



Pmf = Ch[pmf] + CdciPmf] + jOdniPrnf]- (64) 

The Hamiltonian superoperator is given by: Ch[pmf] = —jiHsjPmf], where H = Hg = Hm + Hf + Vmf is a JCM 

type Hamiltonian, Hm = ha (E) i-f, Hf = tiujla ® a^a; Vaf = hX[a2i (8) -|- a^^ a]. CdC'[Pmf] and CdH[Pmf] are 
the dissipative cold and hot Lindblad super operators, respectively: 

^dclPrnf] = ro2{(no2 + l)(['7'02/3Tri,/,Crj2] + [o-02,Prr!,/Crj2])+«02(['7'j2PrTi./,Cro2] + [o'j2jPrn/Cro2])} 

jCdH[Pmf] = ^oi{{mi + i){[(ToiPrnf,crli] + [croi,pmf(Tli])+nQi{[alj^Pmf,(Toi] + [(Tli,Pmfcroi])}, (65) 

where and Fqi are the Weiskopf-Wigner decay constant associated with the cold and hot reservoirs, respectively, 
and nQ2 and uqi are the number of thermal photons in the cold and hot reservoirs, respectively. The temperature of 
the thermal photonic reservoirs is given by: 

The matrix form of the matter operators is given by: 

1 \ / 1 

(721 = I I (701 = I (702 = 

/ \ 
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The extended dissipative JCM master equation (eq. I64p can be obtained by summing the Hamiltonian contribution 
and the two dissi pat ive contributions. Alternatively, it can be derived in a similar fashion to the dissipative JCM 
master equation Assuming that the nature of the interaction between each pair of matter levels is dipole 

coupling, due to parity considerations not all three transitions would be allowed by dipole coupling. However, this 
issue is avoided in systems with a break in symmetry [2l| . 
Let us examine the energy flux of the matter: 

Em = Tr„{p„i7„i} = --Tr{p^/[i?^, Vmf]} + Tr{Cd[pmf]Hm}, (67) 

where Hm is the matter Hamiltonian without the tensor product with 1 f . By substitution into eq. I43i heat flux and 
power for the matter are given by: 

Qm = T^^{i^dc[pmf]Hm} +'^T:{C-dH[pmf]Hm} = QmC + QmH 

Pm = -■^Tr{p^/[Jf^,y^/]}, (68) 
where QmC(H) = ^T:{i^dC(H)[Pmf]Hm}- The energy flux of the selected cavity mode (field): 

Ef = TvfipfHf} = -'-TT{prr,f[Hf,Vrr.f]} = Pf , (69) 

where Hj is the field Hamiltonian without the tensor product with andP/ is the power associated with the field. 
The energy flux of the full matter-field system is given by: 

Emf = Tr{Cdc[pmf]H} + Tr{CdH[p^f]H}, (70) 

where each heat flux component is a sum of two contributions: 

^T^{^dC{H)[Pm.f]H} = Tl {CdCiH)[Pm.f]{Hm + Vmf)} = QmC{H) + QvC{H), (71) 

where Qvc(h) = ^''^{^dC(H)[Pmf\V}- Since we are in perfect matter-field resonance, P,„ = — P/ [[Hm-,Vmf] = 
-[Hf,Vmf]), hence: 

Em + Ef — Emf — Qv- (72) 

Qv vanishes if the off-diagonal matrix elements of Pmf are purely imaginary. Note that to an observer looking on the 
matter alone, work flux (power) and heat fluxes correspond to the traditional view of the first law of thermodynamics 
in which energy is divided into work and heat. The field which is the work source either receives or emits energy to 
the working medium (the matter) in the form of power. 

In a forthcoming publication we give a full dynamical and thermodynamical analysis of the ED- JCM, and show 
that it acts as a quantum optical amplifier [l^ . 

IV. CONCLUSION 

In this paper we have considered the thermodynamics of unipartite systems (systems coupled to an external time 
dependent field and heat reservoirs), as well as bipartite systems. In the latter case, a supplementary part of the 
system replaces the time dependent field. 

For unpartite systems, we gave a minor generalization of Alicki's definitions for heat flux and power, by extending 
these definitions from the Schrodinger to the Heisenberg and interaction pictures. In this generalization, partial time 
derivatives replace full time derivatives in the definitions for heat flux and power. We then presented an alternative 
approach to the partitioning of the energy flux into heat flux and power in which the interaction energy is not included 
directly, that presages our bipartite treatment. We showed that at steady state, if the interaction energy contribution 
to the heat flux vanishes, this alternative partitioning leads to a Carnot formulation of the second law. 

We then turned to bipartite systems. We presented a novel definition for power, based on the energy fluxes of the 
individual subsystems, that is the natural generalization of our alternative partitioning of energy flux in unipartite 
systems. The first law of thermodynamics was derived in differential form in two different ways. The first relies on 
looking at each subsystem independently, while the second relies on looking at the full bipartite system. Although 
any partitioning of the energy fiux is consistent with the first law, the ultimate test of a good partitioning is the 
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fulfillment of the second law. For the partitioning into (full bipartite system) + (reservoirs) the second law in the 
Clausius form follows almost trivially from Spohn's entropy production function. However, for the partitioning into 
(subsystem A) + (subsystem B + reservoirs) we found that there is generally no second law of the Clausius type due 
to oscillations in the partial entropies of the subsystems. Nevertheless, at steady state both the Clausius and Carnot 
versions of the second law are satisfied under some well-controlled conditions. In a forthcoming publication we show, 
both analytically and numerically, that for the ED-JCM both the Clausius and a Carnot versions of the second law 
are satisfied at steady state [171] . 
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